Coexistence of charge density wave and spin-Peierls orders in quarter-filled quasi-one 

dimensional correlated electron systems 
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Charge and spin-Peierls instabilities in quarter-filled (n = 1/2) compounds consisting of coupled 
ladders and/or zig-zag chains are investigated. Hubbard and t-J models including local Holstein 
and/or Peierls couplings to the lattice are studied by numerical techniques. Next nearest neighbor 
hopping and magnetic exchange, and short-range Coulomb interactions are also considered. We 
show that, generically, these systems undergo instabilities towards the formation of Charge Density 
Waves, Bond Order Waves and (generalized) spin-Peierls modulated structures. Moderate electron- 
electron and electron-lattice couplings can lead to a coexistence of these three types of orders. 
In the ladder, a zig-zag pattern is stabilized by the Holstein coupling and the nearest-neighbor 
Coulomb repulsion. In the case of an isolated chain, bond-centered and site-centered 2k f and 4A:f 
modulations are induced by the local Holstein coupling. In addition, we show that, in contrast to 
the ladders, a small charge ordering in the chains, strongly enhances the spin-Peierls instability. 
Our results are applied to the NaV^Os compound (trellis lattice) and various phases with coexisting 
charge disproportionation and spin-Peierls order are proposed and discussed in the context of recent 
experiments. The role of the long-range Coulomb potential is also outlined. 

PACS: 75.10.-b, 75.50.Ee, 71.27. +a, 75.40.Mg 



I. INTRODUCTION 

Quasi-one dimensional correlated electrons systems at 
the commensurate filling of n — 1/2 (quarter-filled band) 
show fascinating physical properties. A widely studied 
class of such materials are the so-called organics charge 
transfer salts like the Bechgaard salts ||. These sys- 
tems which consist of stacks of organics molecules form- 
ing weakly coupled one dimensional chains exhibit, at 
low temperature, a large variety of exotic phases such as 
superconducting, spin density wave, charge density wave 
(CDW) and spin-Peierls (SP) phases. 

The vanadium inorganic compound NaV205 is also be- 
lieved to be a nearly perfect realization of a quarter-filled 
low dimensional system. Therefore, the nature of the SP 
phase below T$p — 35K is expected to be quite dif- 
ferent from the one occurring in the more conventional 
antiferromagnetic Heisenberg chain CuGe03 (Ref. Q). 

The NaV2C>5 system is built from weakly coupled 
planes whose structure is shown in Fig. |l|. It can be 
depicted as an array of parallel ladders (Fig. |l|(b)) cou- 
pled in a trellis lattice. Note that a small buckling of the 
V plane can be neglected in first approximation. Oxy- 
gen atoms (not shown) are located at the center of the 
vertical and horizontal bonds of Fig. |l](a) and lead to ef- 
fective hopping matrix elements and antiferromagnetic 
(AF) super-exchange interactions. LDA band struc- 
ture calculations Q and estimations based on empirical 
rules lead to similar values of the hopping amplitudes 
along and perpendicular to the ladders, in ~ 0.15eV and 
t± ~ 0.35eV respectively. However, some controversy 
remains regarding the magnitude of the diagonal hop- 



ping (see Fig. |l|(c)) t xy with values ranging from 0.012eV 
(Ref. D) to 0.3eV (Ref. §). 

Although the average valence of the vanadium in 
NaV205 is 4.5 (half an electron per vanadium d-orbital 
on average), the exact nature of the charge ordering is 
still under active debate. Early X-rays diffraction experi- 
ments [Q were pointing in favor of a non-centrosymmetric 
structure implying two inequivalent vanadium sites. It 
was further suggested that magnetic spin-1/2 V 4+ 
were forming one-dimensional (ID) chains separated by 
non-magnetic V 5+ chains. Based on this assumption, 
various theoretical analyses of this material were at- 
tempted using dimerized Heisenberg chains J8j or spin- 
phonon models ||. 

Recently, a new structure refinement at room tempera- 
ture |l(]] suggested that, in contrast to earlier reports 0, 
cv'-NaV205 would have a centrosymmetric space group 
implying only one kind of vanadium site. Besides, below 
the transition temperature T$p, joint neutron and X-ray 
diffraction experiments JTT[ | reveal new superlattice re- 
flections which can be ascribed to a lattice modulation 
associated to displacements of predominantly V atoms. 

The insulating character of NaV^Os could, in fact, be 
simply understood in the framework of the quarter-filled 
Hubbard or t-J ladders ||[l2 13 without invoking any 
charge order mechanism (as it is expected at room tem- 
perature). A simple analytic picture, valid when the hop- 
ping along the legs is small |^,^2|, gives a finite charge 
gap of the order of (for a very large Hubbard repulsion U) 
the splitting 2t± between the bonding and antibonding 
states of the rungs. The existence of a metal-insulator 
transition in quarter-filled t-J ladder has been confirmed 
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by numerical calculations of the single particle density of 
state Q and the charge gap Jl3| for realistic parameters. 
In the insulating state, the ground state (GS) configura- 
tion corresponds to the single occupancy of the bonding 
states on each rung and it has been argued that, in this 
case, the low-energy processes can be described by an ef- 
fective AF Heisenberg chain fs|]l5|], hence giving some 
relevance to earlier descriptions of this material H|). 
Furthermore, angle-resolved photoemission data at room 
temperature fl4| are consistent both with a description 
in term of an effective half-filled t- J chain Q or in terms 
of a quarter-filled t- J ladder |L2| . 
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FIG. 1. (a) Structure of the vanadium planes of NaV^Os 
showing ladder (b) and chain (c) patterns. Each vanadium 
carries, on average, a charge of +0.5. 



Recently, motivated by the experimental studies of 
NaV2C>5 at low temperature, charge ordering has been 
investigated theoretically p^ , p7| , ^3| , ^5[ . It has been ar- 
gued that repulsion between electrons on neighboring 
ladders can lead to charge disproportionations where 
the electrons are localized on a single leg of each lad- 
der jl3) (as in the original structure proposed in Ref. R ) 
or form a zig-zag pattern within the ladder |H^ , |i7| . |l5 1 . 
Although the Madelung energy, which includes also a 
long-range (LR) Coulomb interaction, slightly favors the 
chain charge ordering, it has been suggested [|l5| that, 
in contrast, the underlying coupling with the lattice 
could stabilize the zig-zag CDW. Experimental features 
in the optical conductivity at intermediate energies (0.6- 
2.5eV) Q could be reproduced in a calculation of the t-J 
model on the trellis lattice only by assuming an ad- 
hoc charge dispro por tionated GS, although low-energy 
charged magnons [|1 8f| could not be found. 

In this paper, our aim is to investigate further, by nu- 
merical exact diagonalization (ED) techniques, the in- 
terplay between the electron-electron Coulomb repulsion 
and the electron-phonon (or spin-phonon) couplings. In- 
deed, charge disproportionation could come from a Hol- 
stein coupling, possibly related to the apical oxygen, as 
well as from a LR Coulomb repulsion. Although the em- 
phasis is put here on the understanding of the physical 



properties of NaV205, we believe our results can also be 
relevant for organic salts. The likely smallness of the 
interladder as well as many of the experimental results 
point to the possibility that the charge or SP instabilities 
are driven primarily by the ladder physics. |^,^| Thus, 
the first step of our analysis consists of the study of these 
structures by considering isolated quarter-filled t-J lad- 
ders including lattice-local charge couplings (Holstein) 
or modulations of the bond parameters (Peierls and SP 
couplings). We shall then try to determine whether the 
states found for the ladders are also consistent with the 
physics of the chains. Eventually, we shall consider the 
trellis lattice (Fig. |l|(a)) of NaV20s where the most likely 
states will be determined based on the previous results 
obtained for the individual ladder and zig-zag chain sub- 
lattices. 



II. ISOLATED LADDER: CHARGE INSTABILITY 

We shall first consider the quarter-filled anisotropic t-J 
ladder (Fig.@(b)) in the presence of an Holstein-coupling, 

H = H u + H V + H H (1) 
Htj — J|| ^~^(S^q • Sj + i iCK — — ni ja rii+i mCt ) 

i,a> 

+ jj_y^(Sj ? i ■ s 1;2 - -^,1^,2) 

i 

+ i_L^(cj,l;<A2;<r + h.C.) 

Hv = ^^©«M n m,i + ^,2^+1,2 + ^,1^,2) 

i 

h h = Y n ita 5 iia + -K ^ s i,a 1 



where c\ a . a = Ci y0l[ _ a (l — n i a . a ) are hole Guzwiller pro- 
jected creation operators (the large on-site Coulomb in- 
teraction prevents doubly-occupancy) and the index a 
stands for a chain index (— 1,2). We consider a priori 
different fermion hopping amplitudes (tii, t±) or mag- 
netic exchange interactions (Jm, J±) along the legs and 
along the rungs. Hamiltonian (!]) can be viewed as the 
strong coupling limit of a Hubbard ladder so that one 
can assume a relation of the form J\\/J± = (t\\/t±) 2 be- 
tween the parameters. A nearest neighbor (NN) Coulomb 
repulsion V has been included. Note that the electron- 
phonon coupling has been absorbed in the definition of 
the lattice displacement Si. The magnitude of the cou- 
pling to the lattice is then given by a single parame- 
ter namely the inverse of the lattice stiffness 1/K. The 
phonons have been given an infinite mass since the charge 
and spin dynamics are assumed here to involve smaller 
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time scales than lattice fluctuations (adiabatic approx- 
imation). Note that the on-site displacement Si corre- 
sponds in fact to an effective parameter which might 
combine several effects. Physically, it could correspond 
to the displacement of the apical oxygen and/or to the 
displacements of the neighboring in-plane oxygen atoms 
toward or away from the V site. 

Before proceeding with the study of Hamiltonian (Q) , 
it is instructive to recall the properties of the t-J ladder 
in the absence of electron-phonon coupling and NN re- 
pulsion. The electronic properties of this model at quar- 
ter filling have been investigated previously |||l2 13 and 
the existence of a metal-insulator transition has been 
shown |13,0]. It is believed that the system becomes 



insulating for (approximately) t± > 2t\\. Physically, this 
corresponds to the situation where the bonding and anti- 
bonding bands are completely separated, the lower band 
becoming effectively half-filled so that an arbitrary small 
repulsion opens a charge gap. This physical situation 
might be relevant for the insulating phase of the NaV2 0s 
material. 

Although the lattice is considered here in the adiabatic 
approximation, no supercell order is assumed a priori 
and the lowest energy equilibrium lattice configuration is 
obtained through a self-consistent procedure (as in the 
rest of the paper). Indeed, the total energy functional 
E({Si^ a }) can be minimized with respect to the sets of 
distortions {Si^ a } by solving the non- linear set of local 
coupled equations, 



K8 i>a + (rii,a) = 0, 



(2) 



where (...) is the GS mean value obtained by ED (using 
the Lanczos algorithm) of Hamiltonian ([!]). Since the 
second term depends implicitly on the distortion pattern 
{Si^ a }, Eqs. (J2j) can be solved by a regular iterative pro- 
cedure. @ 

The phase diagram of the Holstein-t-J ladder model is 
shown in Fig. ^(a) for a realistic set of parameters. These 
results have been obtained by studying 2x6 and 2x8 
clusters. In the absence of NN repulsion V, a rapid tran- 
sition occurs from a uniform phase (U) at small coupling 
(or equivalently large lattice rigidity) to a localized phase 
at large coupling. This strong coupling phase is charac- 
terized by a charge ordering with two types of (almost) 
completely empty or completely occupied sites arranged 
in some disordered patterns. More interestingly, a CDW 
phase with a "zig-zag" arrangement of the excess charge 
(which exists also for V = only in a very narrow region 
around 1/K ~ 2.5) is stabilized by the NN repulsion V. 
In contrast to the localized phase, the charge dispropor- 
tionation in this zig-zag CDW state is not complete. No- 
tice that there is a finite critical value of 1/K associated 
to the stability of the CDW phase. This feature might 
be due to the fact that the uniform ladder for 1/K = 
is, for an anisotropy ratio of tn/t± = 0.4, already in an 
insulating state with a charge gap (see Ref. fill). 



III. ISOLATED LADDER: COEXISTING 
CHARGE AND SP ORDERS 

The possibility of a coexisting SP order in the previ- 
ous CDW state can be studied by considering additional 
Peierls and SP couplings realized by making the following 
substitutions in Hamiltonian (111) ; 



*n - *n(i + 0> 



Jll 



J\\(l + gS, 



B ) 

1 . ex. / 



(3) 

(4) 



and by adding a new elastic term \Kb J2i a(^Pa) 2 - 

Note that the spin-Peierls order is, for a quarter-filled 
band, intrinsically linked to a bond order wave (BOW) 
characterized by a modulation of the hopping ampli- 
tudes j2lj as given by Eq. (ph. In other words, the ex- 
pectation values (cj Q.g.Cj+i^o-) and (Si )a -Si+i, a ) should 
exhibit similar modulations along the legs. For simplicity, 
we shall assume here that the magnetic exchange intera- 
tion Jy involves two virtual hops tu so that one can take, 
in first approximation, g ~ 1. However, the qualitative 
results of this study do not depend on the choice of g. 
It should be noticed in addition that we shall consider 
modulations only along the legs of the ladders since our 
analysis showed that modulations of the rungs were not 
favored. 

Since our purpose is to investigate the role of the ad- 
ditional Peierls coupling on the zig-zag CDW state, we 
shall first consider an ad-hoc effective potential, 



err 



(xi. c ,+yi, a ) 



(5) 



replacing the terms Hjj + Hy in Hamiltonian ([!]) in or- 
der to stabilize the zig-zag pattern. Qualitatively, this 
term mimics the combined effects of the Holstein lattice 
coupling and the NN Coulomb repulsion. The previous 
method based on Eq. (|J) can be extended here to deal 
with the lattice-bond interaction. The lowest energy con- 
figuration is found by minimizing the GS energy with re- 
spect to the full set of parameters {Sf a }. This is achieved 
by solving exactly, using the same iterative procedure as 
before, the set of non-linear equations, 



^B^fa + <?J||(S, lQ • Si+i )Q — -n^ a ni + i^ a ) 
+ hi^ia^i+x,^ + h.c.) =0. 



(6) 



on finite clusters. As shown in the phase diagram of 
Fig. ||(b), a finite value of the coupling strength 1/Kb ~ 
2 is needed to produce a dimerized pattern overimposed 
to the CDW state. Above this critical value, a finite 
modulation of the form 8f a = 5 B (—l) Xi ' a (the ladder 
is oriented along the x-axis) appears simultaneously with 
the CDW order parameter for any V e g. As can be seen in 
Fig. |||(a), the charge disproportionation, An = — 1|, 
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FIG. 2. Typical phase diagrams of ^--filled anisotropic 
t-J ladders (a) as a function of the NN Coulomb repulsion V 
and the on-site Holstein coupling of strength 1/K; (b) as a 
function of the Peierls lattice coupling and the effective po- 
tential V c fi ; (c) as a function of Holstein and Peierls lattice 
couplings. These results have been obtained by ED of small 
clusters with anisotropy ratios t^/t± = 0.4 and J\\/J± ~ 0.16. 
Shaded regions are unphysical (see text). 

increases almost linearly with the magnitude of the effec- 
tive potential and depends weakly on the magnetoelastic 
coupling as expected. On the other hand, the amplitude 
5b of the dimerization (Fig. ||(c)) increases strongly with 
I/Kb when I/Kb > 2. However, for very strong cou- 
plings, our method gives an unphysical solution 5b > 1 
(implying ferromagnetic bonds) shown by a shaded re- 
gion in Fig. |^(b) (as in the rest of the paper). This sig- 
nals that our simple initial model breaks down for such 
a strong coupling with the lattice and that other terms 



should be included (quartic terms, etc.). Interestingly 
enough, the presence of a zig-zag charge pattern is not at 
all required to observe a dimerization in the ladder. On 
the contrary, as it is clear in Fig. [j](c), the dimerization 
occurs for V c s — (i.e. for a uniform distribution of the 
charge) and is even weakly suppressed by V e s- 
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FIG. 3. CDW order parameter An (a,b) and bond modu- 
lation 1 5f | (c,d) of a i-filled t-J ladder for various lattice bond 
stiffness Kb (as indicated on the plots). (a,c) as a function of 
the effective potential V4a; (b,d) as a function of the on-site 
Holstein coupling strength 1/K. In the "localized" phase (on 
the right of the dashed line), estimates based on an average 
over the system are shown. These results have been obtained 
by ED of small clusters with anisotropy ratios t\\/t± — 0.4 
and J\\/J±_ = 0.16. 

We now turn to the investigation of the full Hamilto- 
nian including both Holstein and Peierls-like lattice cou- 
plings, the amplitudes of each of them being character- 
ized by the two independent coupling constants 1 /K and 
1 /Kb respectively. The GS configuration for the site and 
bond deformations can be obtained exactly on small clus- 
ters by solving simultaneously, by the previous recursive 
method, the two sets of non- linear equations (||) and (||). 
As before, at each step of the recursive procedure, the 
quantum mechanical Hamiltonian is diagonalized by us- 
ing a Lanczos algorithm. Our results, for the same set of 
parameters i M , J,, as above, are summarized in the phase 
diagram of Fig. |2|(c) with the corresponding behaviors of 
the local CDW order parameter and bond deformation 
shown in Figs. ||(b,d). In contrast to the previous model, 
no zig-zag CDW ordering appears at small Holstein cou- 
pling 1/K. A uniform dimerized (D) phase is stable for 
large enough 1/K B - It is important to notice that the 
dimerized zig-zag (D-ZZ) phase is found only in a narrow 
region located around 1/K ~ 2.5 while a larger Holstein 
coupling stabilizes a localized phase (with no well-defined 
periodic pattern). In addition, these data suggest that 
the Peierls coupling rather tends to suppress the charge 
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order and hence to destabilize the D-ZZ phase with re- 
spect to the D phase. However, the region of stability 
of the D-ZZ phase (of particular interest in the case of 
the NaA^Os material) is expected to be greatly enhanced 
by a NN Coulomb repulsion (see Fig. 0(a)) as we have 
checked numerically. 



IV. ISOLATED CHAIN: CHARGE INSTABILITY 

As mentioned above, in NaV2 05, the ladders are cou- 
pled by diagonal bonds forming zig-zag chains. Besides, 
as mentioned in the Introduction, many quasi-lD molec- 
ular compounds contains weakly coupled chains. There- 
fore, we shall now consider the role of the Holstein and 
Peierls lattice couplings on one-dimensional chains. For 
completeness, we start with a ID Hubbard model at 
quarter-filling [n = 1/2) coupled with an on-site (clas- 
sical) phonon field, 

H = t y^(c| !0 .Ct + l; g + h.C.) + U )>"//,: //,:, 

+ v mm+i + nA + \ K J2 5 i ■ ( ? ) 

i i i 

At quarter-filling, the Fermi wave vector is given by 
Q2k F = ■§ so that, at small U , one expects an instability 
towards a 2fcp-CDW state of wavevector X2k F = 4a (a is 
the lattice spacing) mediated by the electron-phonon cou- 
pling. In contrast, for large U, the system becomes more 
similar to a gas of interacting spinless fermions and the 
instability is likely to occur at wavevector 2k F F = AUf- 
More generally, we can parametrize the charge density 
as,© 

1 r 

(rii) = - + A 2kp cos (27T— + $ 2 fe F ) 
I 4a 

+ A 4fcF cos(27r^ + $4fc F ). (8) 
la 

A schematic phase diagram is shown in Fig. ^(a) for V = 
0. The uniform U metallic phase (A 2 k F = A^ F = 0) is 
restricted to a region at small electron-phonon coupling. 
Above a critical line 1/K vs U, three different insulating 
CDW phases can be distinguished; (i) at small U, a 2kF~ 
CDW phase (A^ ~ 0) centered on the sites, i.e. with 
&2k F = 0; (ii) at intermediate U (in the range 4-8), a 
frond-centered 2/c_f-CDW phase i.e. with $2fc F = f ; (hi) 
at large U, a 4fc F -CDW (A 2 k F — 0, $4/^ =0). As seen in 
Fig. |j(b) , a small NN repulsion suppresses completely the 
intermediate phase and enlarges the region of stability of 
the 4fc F -CDW phase |§. 

We now restrict ourselves to the strong electron corre- 
lation limit (t-J model) to discuss further the role of the 
NN repulsion V. Besides, since in NaV205 the zig-zag 
chains form part of the trellis lattice we include NNN 
hopping matrix elements t' and exchange integrals J' 



which correspond to the interactions along the legs of 
the ladders. Hence, the Hamiltonian is, 

H = H u + H V + H H (9) 

H t J = J X^( Sl ' Si +! _ -^ n i n i+l) + t X!( 5 U^+ 1 ' CT + 

i i,cr 

+ J'y^(S, • S i+2 - ^nin l+2 ) + t'^2(cl. t(T Ci + 2;a + h.c.) 

i 
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FIG. 4. Typical phase diagrams of a ^-filled Hubbard 
chain in the presence of an on-site Holstein coupling with- 
out (a) or with a nearest neighbor Coulomb repulsion V (b) 
obtained from ED of small periodic chains. 

The phase diagram of Fig. |5|(a) corresponding to the 
J' = t' = case confirms the stabilization of the AUf 
charge order by the repulsion V . It can be seen that, sim- 
ilarly to what happens in the t-J ladder case (Fig. ||(a)), 
there is a phase of charge disproportionation occurring 
in the region of large Peierls coupling if -1 , with a defi- 
nite pattern. Some important qualitative changes in the 
phase diagram are introduced by the NNN hoppings and 
exchange interactions t' and J' as seen in Fig. ||(b) with 
the appearance of the bond-centered 2/c F -CDW phase 
stable at intermediate Holstein coupling. As in the Hub- 
bard chain, this phase is suppressed by an intermediate 
NN Coulomb repulsion V. 
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FIG. 5. Typical phase diagrams of a ^- -filled t-J chain 
for J = 0.3 as a function of the NN repulsion V and the 
strength of the on-site Holstein coupling 1/K obtained from 
ED of small periodic chains. A hopping t = t and a mag- 
netic exchange J' = J between NNN sites have been included 
(omitted) in (b) (in (a)). 



V. TRELLIS LATTICE: CHARGE INSTABILITY 

At this stage, it is interesting to apply our results 
(valid, strictly speaking, in the case of isolated chains 
or ladders) to investigate the possibility of charge order- 
ing m the 2D trellis lattice of NaV 2 5 (see Fig. 0(a)). 
Obviously, the type of order the most likely to appear 
depends whether the chains or the ladders are the main 
structures of this compound. Since the interladder cou- 
plings, although still not well determined, are quite likely 
to be smaller than the couplings along the rungs of the 
ladders it is to be expected that the physics of the lad- 
ders is going to dominate the behavior of the trellis lat- 
tice. However, this possibility has not been proven so 
far. Thus, here we shall rather try to determine the most 
probable 2D order compatible with both the chain and 
the ladder physics. A 4fcf charge order in the zig-zag 
chains would imply the 2D formation of parallel chains 
of average charge 0.5 + An and 0.5 — An. This corre- 
sponds to the structure originally proposed by Galy and 
coworkers [Q (assuming a complete disproportionation 
An = 0.5). However, experimental evidences accumu- 
late in favor of a modulation of the charge also in the x 
direction of the ladder legs, hence more compatible with 
the 2&i?-CDW orders in the zig-zag chains with NNN cou- 
plings. Two typical patterns exhibiting bond-centered or 



site-centered CDW along the zig-zag chains are shown 
in Figs. |6|(a) andj|(b) respectively. Note that only the 
structure of Fig. n(a) is fully compatible with the zig-zag 
CDW ordering in all the ladders. Based on our previ- 
ous studies, we believe this structure should be stable 
in the 2D trellis lattice for large on-site Coulomb repul- 
sion U, intermediate short-range Coulomb repulsion and 
when the hoppings and magnetic couplings along the legs 
(J|| = J', <|| = t') exceed the cross-bond terms J xy and 
t X y (corresponding to J and t in the chain model). It is 
interesting to notice that, in our picture, the charge or- 
dering is not saturated (An < 0.5). However, a doubling 
of the periodicity occurs in the direction of the ladders for 
any finite value of the order parameter An consistently 
with recent experiments flllfl . 




(a) 



(b) 



FIG. 6. Typical patterns of charge ordering on the 
NaV205 lattice structure showing a supercell of 4 sites along 
the zig-zag chain. The filled (open) symbols correspond to 
an excess (depression) of charge compared to the average 
charge of 1/2 (grey sites have a density close to the aver- 
age density), (a) Bond-centered charge density wave of pe- 
riod 4 (Q = 2/cf) with two types of non- equivalent sites, (b) 
Site-centered charge density wave of period 4 (Q = 2kp) with 
three types of non-equivalent sites. 

So far, only short range Coulomb repulsion has been 
considered. We briefly comment on the role of the long 
range Coulomb repulsion. We have done simple classical 
calculations of the Madelung energy corresponding to the 
various charge configuration introduced before. As previ- 
ously mentioned, most charge ordered states have lower 
energies than the uniform state. Therefore, we believe 
that the long range part of the Coulomb repulsion should 
increase further the tendency towards charge dispropor- 
tionation. 



VI. ISOLATED CHAIN: COEXISTING CHARGE 
AND SP ORDERS 

We now proceed with the investigation of the last issue, 
namely the role of an additional Peierls coupling in the 
chains. We shall follow a method similar to the one used 



G 



before for the ladder structure. For simplicity, we use 
an effective potential to stabilize the bond-centered 2kp- 
CDW phase found previously (see Fig. ||(b)), 



H e g = V e s ^2 Ui C0S ( 27F ^ + f ) • 



(10) 



replacing the terms Hy + Hh of Hamiltonian (^|). This 
effective potential can also include the effects of the LR 
Coulomb interaction on the trellis lattice discussed in the 
previous section. The amplitude of the CDW can then 
be tuned by V e s- We first consider a Peierls coupling in- 
troduced in the NN bonds by the following substitutions, 



t 



J 



- t(l 
J(l H 



2Sf), 



(11) 
(12) 



where the corresponding elastic energy has the usual form 
^KBj2i(^F) 2 ■ GS properties are obtained by solving 
iteratively the set of non-linear equations, 



K B 5f + 2J(Sj • S, ;+ i - 

+ t(cl. a Ci + i ;(7 + h.c.) = . 



(13) 



Typical phase diagrams are shown in Figs. 0(a) and 
^(b) for a simple t-J chain and a t-t'-J-J' chain respec- 
tively. Interestingly enough, a spin-Peierls instability 
(with the same periodicity \2k F = 4a as the underlying 
CDW) occurs for arbitrary small coupling 1/Kb- This 
phase is characterized by three types of bonds in the 4- 
sites unit cell; weak (W), strong (S) or intermediate (I) 
bonds depending on the mag mtudesofof. The sequence 
I-S-W-S in the t' = J' = case corresponds in fact to 
a mixture with attj? component. Note that, in the ab- 
sence of V e e and for 1/ Kb > 1, a simple dimer phase of 
sequence W-S-W-S (pure 4kp component) is stable. As 
expected, the external potential generates an additional 
2fcp component whose magnitude increases with V c g . It 
is interesting to mention here that similar orders occur 
e.g. in quasi-lD organic conductors such as the SP com- 
pound (TMTTF) 2 PF 6 ; while the metallic phase is dimcr- 
ized, a quadrupling of the unit cell is observed in the low 
temperature insulating phase. On the other hand, in the 
case of the zig-zag chain with NNN couplings, an almost 
pure single- Fourier (2k f) CDW is observed with a S-I- 
W-I sequence (Figs. p](b)). It is important here to stress 
that, in contrast to the ladder case studied previously, 
the SP order occurs only in the presence of charge order- 
ing. Furthermore, the SP order parameter follows closely 
the magnitude of the charge disproportionation as it is 
clear from Figs. ||(a) and ||(c). 

For completeness, we should also assume in the zig- 
zag chain the possibility of an additional Peierls coupling 
acting on the bonds connecting NNN sites. This effect is 
then characterized by a new set of bond variables {8f'} 
and a different elastic constant K' B . In the real materials, 
these new parameters are not completely independent 
from the previous ones. However, since the exact 
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FIG. 7. Typical phase diagrams of a ^-filled t-J chain 
for J = 0.3 as a function of the effective potential V e g (see 
text) and the strength of the Peierls bond couplings 1/Kb and 
1/K' B obtained from ED of small periodic chains. A hopping 
t 1 and a magnetic exchange J' between NNN sites have been 
included (resp. omitted) in (b) and (c) (resp. in (a)). In (a) 
and (b) (resp. (c)) the Peierls coupling is included only on 
the bonds connecting NN sites (resp. NNN sites). 

relations between them, which depend on very fine de- 
tails of the chemistry of the material, are quite difficult 
to determine, we shall here examine the role of this new 
coupling separately from the previous one. It will be easy 
then to add the two effects afterwards, at least at a quali- 
tative level. We then consider the following substitution, 

f-fCl + tff), (14) 
J' -» J'(l + 26?) , (15) 

for the parameters associated to the bonds connecting 
the NNN sites i and i + 2. 
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The set of local equations, 
K' B 5f + 2J'(S l • S 



i+2 



H'H+2 



> 



+ t?(c\ Ci+2;<T +h.c.) = 0. 



(16) 



have been solved as before on small 12 or 16 sites peri- 
odic chains and the results are summarized in the phase 
diagram shown on Fig. |7](c). For any arbitrary small SP 
coupling, a dimerization occurs in the two parallel chains 
formed by the NNN bonds, coexisting with the under- 
lying CDW state. As seen in Figs. ||(b) and ||(d), two 
qualitatively different regimes have to be distinguished; 
(i) for small SP couplings (let say l/K' B < 2), the dimer- 
ization and the charge order appear simultaneously and 
follow a linear behavior as a function of the magnitude of 
the external potential; (ii) for larger couplings, the mixed 
CDW-SP state is stable even in the absence of V c s. In 
that case, the charge disproportionation Art is further 
increased by V c g while the dimerization S B remains al- 
most constant. It is important to stress here that, in 
both regimes, both order parameters follow a very sim- 
ilar behavior as a function of the coupling constants of 
the model suggesting that the two types of orders are 
intrinsically linked. 
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FIG. 8. CDW order parameter An [in (a) and (b)] and 
bond modulation |<5j B | [in (c) and (d)] of a i-filled t-J chain for 
J = 0.3 as a function of the effective potential V^t (see text) 
and for various lattice stiffness Kb (as indicated on the plots) 
obtained from ED of small periodic chains. A hopping t' and a 
magnetic exchange J' between NNN sites have been included. 
In (a) and (c) (resp. (b) and (d)) the Peierls coupling is 
included on the bonds connecting NN sites (resp. NNN sites). 



VII. TRELLIS LATTICE: COEXISTING CHARGE 
AND SP ORDERS 

Based on our previous studies of the coexistence of 
CDW and SP orders in isolated ladders and zig-zag 



chains, we shall now attempt to construct a unified pic- 
ture for the 2D trellis lattice of NaV20s. Although this 
approach is not rigorous, we believe it can provide reliable 
qualitative predictions of the low temperature structure 
of this material. Our starting point is the charge ordered 
pattern of Fig. [| likely to be stabilized by a local on- 
site Holstein coupling and/or a LR Coulomb repulsion. 
As seen before, this structure combines the zig-zag CDW 
pattern of the ladders and the bond-centered 2fcp-CDW 
of the chains. If the SP instability is driven primarily 
by a Peierls coupling in the diagonal bonds t xy _= t and 
J xy = J, one should expect the pattern of Fig. g(a) with 
the I-S-W-S sequence along the zig-zag chains. However, 
as argued previously, the stability of the bond-centered 
2&F-CDW in the chains which we assume here requires 
rather sizable NNN couplings t' and J' corresponding to 
the tn and Jn parameters of the trellis lattice. There- 
fore, we believe that the Peierls coupling corresponding 
to these bonds can no longer be neglected. In this case, 
one can argue that the structure shown in Fig. ^|(b) is 
more stable than the one shown in Fig. ||(a). It is ob- 
tained by a superposition of the two types of SP orders 
of the t-t'-J-J' chains shown in Figs, f7|(b) and 0(c), each 
of them being mediated by Peierls couplings in the diag- 
onal bonds and the legs respectively. This state is char- 
acterized by a dimerization in the legs and a modulation 
S-I-W-I in the t-t'-J-J' chain. Note that the SP order 
on the diagonal bonds of the t-t'-J-J' chain is somehow 
"locked" to the charge order i.e. the weak (W) bonds are 
always the ones connecting two sites with an excess of 
charge and so on. However, the actual pattern of the 





(a) 



(b) 



FIG. 9. Typical patterns of coexisting lattice distortion 
and charge ordering on the NaV^Os lattice structure show- 
ing a periodicity of 4 sites along the t-t'-J-J' chain. The filled 
(open) symbols correspond to an excess (depression) of charge 
compared to the average charge of 1/2. These patterns show a 
modulation of the bond exchange couplings along the zig-zag 
chain with three types of bonds and two types of sites (2k f 
order of the charge density). Thick, thin and dotted lines 
correspond to strong, intermediate and weak bonds respec- 
tively, (a) Mixed 2kp-4kF lattice distortion, (b) 2&f lattice 
distortion. 
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chains depends on the relative phase between the two 
SP modulations. Fig. ||(b) corresponds to the simplest 
and most natural choice with only two different kinds of 
triangles. Besides, we believe that this possibility has the 
lowest energy. 



A A' A B' A A' B B' 




(a) (b) 

FIG. 10. Same distortion and charge ordering pattern 
along the zig-zag chains as on Fig. ^(b) but with a doubling 
of the periodicity in the transverse direction. A quadrupling 
(doubling) of the 2D unit cell occurs in (a) ((b)). 

The last step in our procedure consists now to refine 
the final structure using our knowledge of the SP insta- 
bilities in an isolated ladder. In fact, by examining the 
ladders of Fig. §(b), we realize that the type of dimeriza- 
tion occurring in the legs do not correspond to any type 
of SP order found before in the isolated ladder. However, 
although still imposing a fixed type of pattern in the legs 
and zig-zag chains, other arrangements are possible in 
the transverse direction leading to qualitatively different 
structures. Since the charge ordering doubles the period- 
icity in the ladder direction, there are two different config- 
urations A and B of the ladder legs, which transform into 
each other by a translation of a unit vector. In addition, 
due to the coexisting SP-BOW order, a reflection with 
respect to a site perpendicularly to the leg leads to two 
non-equivalent configurations A and B. As can be seen 
straightforwardly, the characteristic pattern of the zig- 
zag chains of Fig. |9|(b), can appear in four different ways 
(all related to each other by a lattice translation and/or a 
reflection with respect to a site), namely AA', AB', BA' 
and BB' (the prime refers to the second leg) . Various or- 
ders can be realized by arbitrary sequences in the trans- 
verse direction involving different configurations in each 
primitive cell. Note that each primitive cell contains two 
zig-zag chains but the pattern in the second zig-zag chain 
is somehow imposed its two neighboring ones. The pre- 
vious arrangement shown in Fig. |t](b) corresponds then 
to a periodic arrangement in the transverse direction like 
-AA'-AA'- (equivalent to e.g. -BB'-BB'- and so on). In 
contrast, the ordering patterns of Figs. |l0|( a) and[lO|(b) of 
the form -AA'-AB'- (equivalent to -AA'-BA'-) and -AA'- 



BB'- respectively double the periodicity in the transverse 
direction as well. It is important to stress that only the 
pattern shown in Fig. |l^(b) show the expected mixed 
CDW-SP ordering of the D-ZZ type in all the ladders 
(see phase diagrams in Figs, ^(b) and (c)) . We can then 
argue that the ordering shown in Fig. [10](b) can be sta- 
bilized both by the couplings in the chains and in the 
ladders. This state corresponds in fact to a check-board 
pattern with a doubling of the unit cell and is consistent 
with recent light and neutron scattering experiments [ pT| . 

VIII. CONCLUSIONS 

To summarize, the role of Holstein and Peierls electron- 
phonon couplings as well as magneto-phonon (SP) cou- 
plings has been investigated in the adiabatic approxima- 
tion in quarter-filled one-dimensional chain and ladder 
systems. A numerical method based on ED techniques 
supplemented by a self-consistent procedure has been 
used to determine various phase diagrams as a function 
of the strengths of the lattice couplings. We have shown 
that, generically, CDW, BOW and spin-Peierls orders co- 
exist in such systems. Moreover, in many cases, the SP 
instability is enhanced by a charge disproportionation. 
Eventually, we have considered the case of the trellis lat- 
tice of NaV2C>5. Based on the previous results for isolated 
chains and ladders, we have proposed a mixed CDW- 
BOW-SP ground state with a zig-zag charge pattern in 
the ladders and a modulated bond structure involving 7 
different bonds. Such ordered state doubles the unit cell 
introducing a check-board type pattern. 

Whether the charge ordering and the BOW-SP tran- 
sition occur simultaneously at the same temperature is 
not yet clear but could be resolved experimentally by 
thermodynamical measurements. In addition, such tran- 
sitions should have different spectroscopic signatures due 
to profound changes in the GS excitation spectrum. On 
one hand, charge ordering is expected to strongly en- 
hance the charge gap (although the insulating character 
of the material at high temperature could be accounted 
for by a local on-site repulsion alone). On the other hand, 
the opening of a spin gap should be associated to the SP- 
BOW transition. 

Very recent X-ray diffuse scattering studies of 
NaV2 05 [^3| suggest that 3D structural fluctuations are 
important above the transition and, hence, that trans- 
verse interactions (other than magnetic) should be im- 
portant to stabilize the ordered phase. In our description, 
the SP mechanism is intrinsically linked to the occur- 
rence of charge ordering, possibly induced by a coupling 
with the 3D lattice (treated here in the adiabatic ap- 
proximation). Therefore, we believe that the mechanism 
proposed in the present work is, to some degree, weakly 
dependent on the actual magnetic couplings of the trellis 
lattice, which, so far, are not precisely known. 
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